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Abstract 


Application of a quasi-3D approach to the 
aerodynamic analysis of several radial turbine 
configurations is described. The objective was to 
improve the rotor aerodynamic characteristics by 
hub-shroud contouring. The approach relies on 
available 2D inviscid methods coupled with boundary 
layer analysis to calculate profile, mixing, and 
endwall losses. Windage, tip clearance, incidence, 
and secondary flow losses are estimated from cor- 
relations. To eliminate separation along the hub 
and blade suction surfaces of a baseline rotor, 
the analysis was also applied to three alternate 
hub-shroud geometries. Emphasis was on eliminating 
an inducer velocity overshoot as well as increasing 
hub velocities. While separation was never elimi- 
nated, the extent of the separated area was pro- 
gressively reduced. Results are presented in 
terms of mid-channel and blade surface velocities; 
kinetic energy loss coefficients; and efficiency. 
Geometries illustrated are not an exhaustive 
attempt at design optimization. The calculation 
demonstrates a first step for a systematic approach 
to radial turbine design that can be used to iden- 
tify and control aerodynamic characteristics that 
ultimately determine heat transfer and component 
life. Experimentation will be required to assess 
the extent to which flow and boundary layer be- 
havior were predicted correctly. 

Nomencl ature 


e kinetic energy loss coefficient 
P static pressure, N/m 

P‘ absolute total pressure, N/m 

P" relative total pressure, N/m 

r radial distance, cm 

z axial distance, cm 

t ) efficiency 

Subscr ipts : 

0 stator inlet 

1 stator exit 

2 rotor exit (trailing edge) 

3 rotor exit (mixed-out plane) 

id ideal 

t-t total-to-total 

Int ro duction 

The radial turbine continues to be a candidate 
component for small gas turbine engines because of 
its high stage work and potential efficiency ad- 
vantage over the axial turbine. Past efforts aimed 
at describing the characteristics of radial tur- 
bines include the work by Balje* who described 


their performance in terms of Mach number, Reynolds 
number, specific diameter, and specific speed. 
Benson^ followed a similar approach but refined 
some of the geometry characteristics to include 
inlet and outlet blade height and nozzle and rotor 
exit angles. Watanabe 3 experimentally examined 
effects of the vaneless space gap, the area ratio 
between impeller inlet and exit, rotor inlet to 
exit diameter ratio, and clearance. Mizumachi 4 
experimentally studied the effect of blade number 
and exducer blade angle. From measurements of 
flow angle and velocity coefficient, he inferred 
what was happening inside the rotor passage regard- 
ing secondary flows and separation. In general, 
these and other efforts in radial turbine perform- 
ance characterization have been limited to overall 
flow conditions into and out of the turbine and 
have not addressed the flow inside the blading. 

For many practical applications, this is perfectly 
adequate. 

For high-temperature, highly-stressed appli- 
cations, however, the radial turbine's acceptance 
has been delayed because of the difficulties 
associated with fabricating a cooled rotor that 
satisfies life requirements. A significant part 
of these difficulties can be traced to a lack of 
basic understanding of the boundary layer behavior 
through the rotor passage. The lack of a syste- 
matic approach for even the aerodynamic calcula- 
tions aggravates an already difficult heat transfer 
and cooling design problem. The resultant uncer- 
tainty in rotor aerothermodynamics has made radial 
turbine design as much an art as it is a science, 
with each designer relying heavily on some unique 
past experience. 

This paper describes the application and re- 
sults of a systematic quasi-3D approach to the 
aerodynamic analysis of several radial turbine 
configurations. The objective was to improve the 
aerodynamic characteristics of the rotor by hub- 
shroud contouring. Overall losses were estimated 
in all cases, but since no experimental data were 
available, a comparison with data could not be 
made. The approach was to use the coupled MERIDL, 
TSONIC, and BLAYER computer code described in 
Ref. 5. These three codes have been coupled and 
linked with a loss analysis model which has been 
shown to agree well with experimental data for 
axial turbines. In the present effort, the same 
analysis, with some minor modifications, was 
applied to four radial rotor configurations. The 
baseline configuration is an early version of a 
20 cm cooled radial research turbine designed for 
an advanced rotorcraft application of 1607 K inlet 
temperature and a work output of 465 J/g. Because 
of a large predicted separation on the hub and 
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suction surfaces of the baseline rotor, an attempt 
was made to minimize or eliminate the diffusion by 
hub-shroud contouring. While holding the stator 
geometry unchanged, three alternate rotors with 
varying hub-shroud contours were analyzed. Blade 
number and blade geometry were held fixed. Results 
are presented in terms of mid-channel velocity 
profiles, blade surface loadings, extent of sepa- 
ration, efficiency, and kinetic energy loss 
coefficients. 

Method/Approach 

The analysis approach used in this paper is 
essentially the one presented by Boyle, Haas, and 
Katsanis in Ref. 5. The analysis procedure is 
applied to both stator and rotor and is iteratively 
coupled to the calculated losses. First, the in- 
viscid two-dimensional computer program MERIDL is 
used to calculate the flow velocities on the hub- 
to-tip midchannel stream surface with an assumed 
pressure drop due to losses. The resultant stream- 
sheet thicknesses are subsequently used in the 
TSONIC program to obtain solutions on five blade- 
to-blade stream surfaces from hub-to-tip. The 
pressure and suction surface exit static pressures 
are made equal for each TSONIC solution by adjust- 
ing the downstream whirl distribution for the 
MERIDL program. This iteration is repeated until 
the static pressures are equal within some toler- 
ance limit. The resulting quasi-3D solution of 
the flow through the passage serves as input for 
the integral method boundary layer code BLAYER to 
compute the boundary layer growth along pressure, 
suction, and endwall surfaces. Profile, mixing, 
and endwall losses are determined from boundary 
layer parameters using Stewart's analysis. 

Empirical models are used to account for incidence, 
disc cavity, rotor tip clearance, and secondary 
flow losses. A mass-averaged overall loss is cal- 
culated and an equivalent pressure drop is deter- 
mined. The newly calculated pressure drop replaces 
the initially guessed value in MERIDL and the en- 
tire procedure is repeated until the two values 
agree within some tolerance limit. 

Except for some minor modifications to account 
for the radial geometry, the empirical loss models 
used were generally the same as those described in 
Ref. 5. The few exceptions can be summarized as 
follows. The same secondary loss correlations were 
used, but with the original authors' coefficients^’^ 
in the secondary loss equation. Since endwall 
losses are computed independently, only the leading 
edge vortex term was retained. Since the tip 
clearance model used in Ref. 5 accounted only for 
radial clearance in an axial turbine, the results 
from Ref. 9 were used to obtain an effective clear- 
ance that accounted for both axial and radial 
clearances found in a radial turbine. An optimum 
incidence angle for the radial rotor was determined 
from the slip factor correlations developed by 
Wiesner in Ref. 10. The incidence loss was then 
calculated from the deviation angle as in Ref. 5. 

The disc cavity loss model of Ref. 5 was used but 
with a reduced rim radius to simulate a scalloped 
backface for the radial rotor. 

Some notes are also in order regarding the 
boundary layer analysis. The BLAYER code has been 
revised to include the effect of radius change. 

For the present analysis, the boundary layer was 
assumed to be entirely turbulent and no laminar 


solution was calculated. Small initial displace- 
ment and momentum thicknesses were specified. The 
separation criteria was skin friction coefficient 
becoming zero or negative. In cases where separa- 
tion did occur, the smoothing feature in BLAYER 
was first used to locally smooth the free-stream 
velocity solutions obtained from TSONIC. If the 
separation persisted after a specified number of 
smoothings, a BLAYER solution would be generated 
for a revised velocity distribution that elimi- 
nated surface diffusion entirely. As applied to 
the blade suction surface at the hub, the revised 
velocity distribution would retain the strong 
acceleration found near the trailing edge and would 
simply maintain a constant minimum velocity over 
the forward portion of the blade. The inherent 
assumption is that this revised solution is indica- 
tive of a boundary layer thickness at the trailing 
edge if the flow re-attached after an initial sepa- 
ration. Since there is no penalty applied to cases 
that separate and re-attach, it is probably not 
very meaningful to compare such cases on calculated 
efficiency alone. In such cases, however, the 
analysis still provides guidance as to the presence 
and extent of the separation. 

Results 

Baseline Rotor Configuration 

Design requirements . - The turbine design 
requirements are summarized in Table I. Size, 
work, and inlet conditions are representative of a 
high-pressure gas-generator turbine for an advanced 
rotorcraft application. The meridional flowpath 
of the complete stage is shown in Fig. 1(a). There 
are 15 stator vanes and 14 rotor blades. The 
stator exit flow angle was approximately 75°. The 
inducer was not swept and the rotor blading had 
structurally radial blade elements. The blading 
profiles were intended to accomodate cooling and 
reflect thickness distributions based on stress 
considerations. A center bore through the rotor 
hub was included to allow passage of an output 
power shaft. This design feature impacted the 
bore stresses considerably and dictated a maximum 
radius of about 5.6 cm for the solid material por- 
tion of the hub. 

Mid-channel v eloc i ties . - Figure 1(b) shows 
the MERIDL orthogonal mesh used for the baseline 
geometry flow analysis. Mid-channel velocities 
calculated along five streamlines are shown in 
Fig. 2. In the inducer entrance region, the flow 
rapidly accelerates as the area decreases with 
radius in the streamwise direction. As the hub and 
shroud contours begin to diverge and also turn 
axially, the streamwise area variation reverses 
and the flow experiences a strong deceleration. 

This deceleration is most pronounced along the hub 
streamline. In the exducer region, where the blade 
turning decreases the relative throughflow area, 
the flow is strongly re-accelerated. 

Blade surface velocities . - The resulting 
blade surface velocities along hub, mean, and tip 
streamlines are shown in Fig. 3 for the baseline 
rotor. As anticipated, the mid-channel diffusions 
are further accentuated on the blade suction sur- 
faces, particularly in the hub region. 

Boundary layer . - A BLAYER analysis of the 
baseline rotor blade pressure and suction veloci- 
ties was done for five hub-to-tip streamlines. 
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Using the hub TSONIC solution, a BLAYER analysis 
was also done for five streamlines along the hub 
endwall. Figure 4 shows the blade suction and hub 
endwall surfaces in nondimensional form as stream 
function and streamwise distance. The shaded area 
represents the approximate extent of the separated 
region for the baseline configuration. Separation 
was predicted by BLAYER when skin friction coeffi- 
cient went to zero or negative. For all the cases 
along streamlines that separated, revised solu- 
tions, as explained earlier, were obtained. The 
re-attachment point was taken to be where the mo- 
mentum thickness began to decrease in the revised 
solution. As indicated in the figure, flow sepa- 
ration for the baseline rotor occurred on the blade 
suction surface from the hub up to 25 percent 
span. Separation also occurred on the hub endwall 
from the blade suction surface to the 75 percent 
streaml ine. 

Contoured Hub-Shroud Configurations 

Geometry . - In an attempt to minimize or 
eliminate the relatively large separation region in 
the baseline rotor, three alternate rotor configu- 
rations were similarly analyzed. In each case, 
only the hub and shroud contours were varied and 
the same input blade sections were used to gener- 
ate the actual blading within the limits of those 
contours. The stator was left unchanged. The 
emphasis in each case was on eliminating or mini- 
mizing the initial overacceleration of the mid- 
channel velocity in the inducer region as well as 
increasing the minimum velocity along the hub. The 
5.6 cm maximum radius for the solid hub portion was 
kept as an anchor for all the alternate hub con- 
tours. Figures 5 to 7 show the three alternate 
hub-shroud geometries. Part (a) of each figure 
shows the alternate contour compared to the base- 
line and part (b) shows the resulting MERIDL mesh. 
The contour A geometry in Fig. 5 has two modifica- 
tions to the shroud contour. To compensate for the 
rapid flow acceleration in the inducer region, the 
flow area in that region was increased by modifying 
the shroud contour. To increase the minimum flow 
velocities along the hub, the entire exit flow area 
was reduced by decreasing the exducer tip diameter. 
In addition to these shroud changes, the contour B 
geometry shown in Fig. 6 incorporates hub contour 
changes to further increase and decrease the flow 
areas in the inducer and exducer regions, respec- 
tively. The contour C geometry shown in Fig. 7 
has one additional modification to the exducer 
shroud contour to further increase the minimum 
velocity along the hub. 

Mid-channel velocities . - The impact of the 
various hub-shroud contour changes on the mid- 
channel velocities can be seen in Figs. 8 to 10. 

The inducer shroud changes of contour A (fig. 8) 
decreased the initial velocity overshoot and the 
additional contouring of the inducer hub in contour 
B (fig. 9) all but eliminated it. The improved 
reaction with reduced exit area raised the minimum 
velocities along the hub for all the alternate 
configurations. The exducer hub contour of rotor 
B (fig. 9) and the exducer tip contour on configu- 
ration C (fig. 10), each increased the minimum hub 
velocities further. It should be noted that the 
variations in mid-channel exit velocities are not 
simply due to continuity responding to area changes 
alone. Closure of the blade loading at the trail- 
ing edge was achieved at five radial locations for 
each rotor by varying the radial work distribution. 


Stage pressure ratios, therefore, were not identi- 
cal and the converged solutions for rotor contours 
A, B, and C had work levels 5.4, 5.3, and 7.9 
percent higher than the baseline rotor, respec- 
tively (see bottom of Table II). Another level of 
iteration involving rotor blade geometry variation 
would have been necessary to achieve identical work 
levels. 

Blade surface velocities . - The corresponding 
blade surface velocities for the contoured rotors 
at the hub, mean, and tip sections are shown in 
Figs. 11 to 13. The trends that were observed in 
mid-channel velocities are similarly reflected in 
hub suction surface diffusion. At some expense in 
increased tip region velocities, contouring the 
hub-shroud did progressively reduce hub suction 
surface diffusion. 

Boundary layer . - A BLAYER analysis as 
described previously was done for the three alter- 
nate rotor hub-shroud contours. Figure 14 shows 
the resulting separated region for contours A, B, 
and C superimposed on the baseline results from 
Fig. 4. The extent of the separation can be seen 
to be progressively reduced. While contour C had 
the smallest region of separated flow, it was never 
entirely eliminated. Although there was some in- 
creased diffusion in the tip region, especially 
for contour C, none of the alternate contours 
showed separation in the tip region. 

Estimated losses . - For completeness, Table 
II presents the loss analysis results. Losses are 
tabulated in terms of kinetic energy loss coeffi- 
cients and total-to-total efficiency. The bottom 
of the table also includes the specific work and 
some pertinent pressure ratios for each rotor con- 
figuration. The stage efficiency numbers represent 
uncooled values. A constant-area exhaust duct was 
included in the analysis because the baseline 
design had been done for a mixed-out downstream 
plane. The stator loss coefficients are constant 
for each configuration, but, because stage pressure 
ratio varied slightly, the efficiency decrement 
due to the stator varied also. Although the over- 
all rotor loss coefficients progressively decreased 
from baseline through contour C, the accompanying 
increasing pressure ratios produced negligible 
differences in overall efficiency. The results in 
Table II are for the average of the two reference 
secondary loss correlations. Using Dunham's cor- 
relation alone resulted in about one-half point 
lower efficiencies. 

Conclusions 

Flub-shroud contouring did reduce the estimated 
extent of the flow separation in a radial rotor. 

The particular hub-shroud geometries illustrated 
do not represent an exhaustive attempt to optimize 
a radial turbine design. Rather, the calculation 
demonstrates a first step of a systematic approach 
to radial turbine design using a quasi-3D approach. 
The significance is that the analysis relies mini- 
mally on past experience and can be used to iden- 
tify and control the aerodynamic characteristics 
that ultimately determine heat transfer and affect 
component life. Experimentation will be required 
to assess the extent to which the analysis predicts 
the correct flow aerodynamics and boundary layer 
behavior. 
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TABLE I. - TURBINE DESIGN CONDITIONS 


Inlet total temperature, K 1607 

Inlet total pressure, N/m^ 1.637x10° 

Specific work, J/g 464.9 

Work factor 1.0 

Rotor tip speed, m/s 682.8 

Rotor tip diameter, m 0.2038 

Rotative speed, rpm 64 000 

Mass flow rate, kg/sec 2.370 

Power, kW 1102 




Stator losses 

Profile and mixing 
Endwall friction 
Vaneless space 
Secondary 
Incidence 

Total 


Rotor losses 

Profile and mixing 
Hub endwall friction 
Secondary 
Incidence 
Tip clearance 
Disk friction 

Total 

Exhaust duct 


Overall uncooled stage eff. 


V p i 

C-d /P 3 

p >; 


Specific work, J/g 


^t-t 


TABLE II. - ESTIMATED LOSSES 


Baseline Contour A Contour B 


e 


An t _ t e An t _ t e &n t _ t 


0.0187 Same Same 

.0165 as as 

•0127 baseline baseline 

.0227 
CL 

0.0706 0.0332 0.0309 0.0310 



J Contour C 

e 

4 Vt 

Same 

as 

basel i ne 

0.0299 

0.0524 

.0031 

.0206 

.0346 

.0875 

.0016 

0.0216 

.0098 

.0248 

.0004 

0.1998 

0.0566 


0.0006 

0 . 

9129 


1.881 1.882 1.882 1.882 

1.324 1.455 1.451 1.535 


3.49 

458.5 


3.78 

483.3 


3.77 

482.7 


3.94 

494.9 











































Figure 2. - Relative mid-channel velocities for baseline 
rotor. 
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(a) Hub (0% streamline). 
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(c) Tip (100% streamline). 

Figure 3. - Relative blade surface velxities for base' 
line rotor. 
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Figure 4. - Extent of separated flow region for baseline rotor. 
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(a) Meridional flowpath. 



(b) MERIDL mesh. 

Figure 5. - Contour A geometry. 
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(b) MERIDL mesh. 

Figure 6. - Contour B geometry. 
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Figure 7. - Contour C geometry. 
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Figure 8. - Relative mid-channel velocities for contour A. 



Figure 9. - Relative mid-channel velocities for contour B. 
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Figure 11. - Relative blade surface velocities for 
contour A. 
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Figure 14. - Extent of separated flow region for baseline and contoured hub-shroud rotors. 
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